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Realistic calculations of nuclear disappearance lifetimes induced by nn oscillations are reported 
for oxygen and iron, using n nuclear potentials derived from a recent comprehensive analysis of p 
atomic X-ray and radiochemical data. A lower limit T n n > 3.3 x 10 8 s on the nn oscillation time 
is derived from the Super-Kamiokande I new lower limit Td(O) > 1.77 x 10 32 yr on the neutron 

© : lifetime in oxysen - Antineutron scatterins lengths in carbon and nicke1 ' needed in trap experiments 

using ultracold neutrons, are calculated from updated N optical potentials at threshold, with results 
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I. INTRODUCTION AND OVERVIEW 



The stability of nuclei, as determined by looking for proton decay P, sets a lower limit also on the lifetime 

' of other processes such as neutron-antineutron (nn) oscillations in free space. This AB = 2 baryon-violating nn 
' oscillation process was pointed out long ago in the pioneering papers by Kuzmin Q , by Glashow Q and by Mohaptra 
and Marshak [6fl. Several quantitative calculations relating the nuclear disappearance lifetime T c i to the nn oscillation 

. n n q 

time T n n have been reported [7|, |8|, |9(. In these calculations a point-like nn coupling 8m = h,T n ^ is assumed. In free 



space, 8m splits the n — n degenerate mass m into mass eigenstates m ± 8m. The nn oscillations between these two 
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00 

mass eigenstates are suppressed in nuclear matter, giving place to decay of neutrons in a stable nucleus. Instead 
of two mass eigenstates one encounters two distinct widths, one which is the nuclear n annihilation width of order 
T n rj 320 MeV for central nuclear densities [10], the other one associated with the lifetime of a bound neutron: 

4 8m 

T d « £p) 8m , (1) 

t n 

where Td = KT^ 1 is the nuclear disappearance width per neutron. These statements follow, more rigorously, from a 
discussion of the temporal evolution of nuclear disappearance owing to nn oscillations [11| . Eq. ([T]) demonstrates that 
nn oscillations are suppressed in nuclei by 31 orders of magnitude, which is the ratio between the time scales 10~ 23 s 
for the strong-interaction n annihilation time r n — KT^ 1 and 10 8 s presumed below for the free-space oscillation time 
T nn = hiSm)^ 1 . In addition to n annihilation, neutrons and antineutrons also experience different nuclear potentials 
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U n and Un, leading within a closure approximation 12| to a refinement of Eq. ([1}: 



^closure ^ (5m) 

< W n > 2 + < (Z7 fl - U n ) > 2 



^closure ^ p rt\ 



where < >= r ft /2 is some appropriately chosen average of (minus) the imaginary part of the n nuclear potential. 
To leading order, from Eq. {T]), the relationship between T nfl and Td is given roughly by 

T nn « 2 (ft Td/T^) 1 ' 2 . (3) 

The values reported for Td in the literature, by convention, are normalized per neutron and take into account secondary 
and absorption processes. Thus, Td essentially stands for the lifetime of a neutron in a stable nucleus. Using the 
best value T d (Fe) > 7.2 x 10 31 yr published by the Soudan 2 Collaboration 3], Eq. Q gives a lower limit estimate 
of T n n > 1.4 x 10 8 s. This estimate is to be compared with the lower limit T nfl > 1.3 x 10 8 s stated by the Soudan 
2 Collaboration using the calculations made long ago by Dover-Gal- Richard [8(. However, this apparent agreement 
might be fortuitous. We comment that the lower limit on T nn derived from nuclear disappearance considerations is 
higher than that determined using nuclear-reactor neutrons directly in searches for nn oscillations. The lower limit 



given by the Grenoble reactor experiment 



3 is T„ 



> 0.86 x 10 8 s. 



Recently, the Super-Kamiokande (SK) collaboration released an improved value for Td in oxygen, Td(0) > 1.77 x 



141 ]. In the present work we report an accurate calculation of the lower limit on r, ift implied by this value of 



10^ yr 

Td(0), using the latest detailed analysis by our group which derived antinuclcon nuclear potentials from antiprotonic 



atom data and radiochemical data 10( . These potentials are essentially isoscalar and apply to antineutrons as well as 
to antiprotons. For completeness, in this study we also calculate the reduced lifetime (see below) in iron which will 
provide the necessary link between a future measurement of X^Fe) and a derivation of a lower limit on T nfl . 

Another experimental method of looking for nn oscillations is to use ultracold neutrons in a trap, first suggested 
by Chctyrkin et al. 



12j | and discussed since then by several other groups, e.g. in Refs. [15|, [16|. A necessary input 
for interpreting such experiments are the nuclear n scattering lengths for ultracold antineutrons in matter. We report 
here on calculations of n scattering lengths for several chosen materials, again based on the analysis of antiprotonic 
atoms. The calculated values of these scattering lengths turn out to be largely model independent owing to the strong 
absorption of low energy antinucleons. 

II. METHODOLOGY 

A. Wave equations, widths and lifetime 

A point-like coupling 8m = hr~^, representing nn, oscillations in free space, connects each of the stationary bound 
single-particle (sp) neutron states to the corresponding n stationary sp state. The sp energies assume complex values 
E v tj = —B^ij — ir„£j/2, where the imaginary part of the energy in the sp state labeled v£j gives the disappearance 
width T^ij for this state. The coupled n-h radial wave equations for the neutron sp wavefunction u„£j(r) and the 
antineutron sp wavefunction to^ (r) are given by 



(r) | 



u v£j( r ) + — 2 — u vt d {r) - U„{r) u v (j(r) - E utj u vtLj {r) + 5m w ut j(r) = , (4) 



_ 2/1 Wul ^ H — Wvl 'j( r ) ~ I c/ n( r ) + i WftW] w fij{ r ) - E vij Wuij{r) + 5m u U £j(r) = , (5) 

where —U n (r) and — (Un(r) + i Wji(r)) are the nuclear potentials exerted by the nuclear core on the neutrons and 
antincutrons, respectively, and fi is the reduced mass. The radial wavcfunctions u v tj(r) and w v ij(r) are regular at 
the origin and decay with r outside of the nucleus. A useful expression for the width T u gj is obtained by multiplying 
Eq. (|4]) by u* £j -(r), and the complex conjugate of Eq. Q by u v £j(r), subtracting the resulting expressions from each 
other and integrating from zero to infinity. The result is 

2 5m Jim {w vtj (r) u* vlj (r)) dr 

Ulj = Fl / \ 1 2 t • W 

The initial-time boundary condition of no antincutrons implies that | w/u \= 0(5m/B), where the binding energy B 
represents any of the strong- interaction entities, such as T n It follows then from Eq. (J6j) that the width T is of 
order (5m) 2 /Tn, in agreement with Eq. fl}. The terms with T and 5m in Eq. (J3]) are negligible, of order (5m/ B) 2 with 
respect to the rest of the terms which coincide with those constituting a stable bound-neutron radial wave equation 

(n) 

in which B v ij = B ul - stands for the neutron sp binding energy. Thus, the solutions u^gj are essentially real functions. 
Eq. ([6]) expresses a relationship between these two minute terms which are neglected below. Turning to Eq. ([5]), all 
terms there are of the same order, except for the T term which is of order (5m/ B) 2 with respect to the other terms 
and, hence, may be dropped off. 

Dropping off the terms of order (5m/ B) 2 in Eqs. (j4} and ([5|), we obtain the radial equations satisfied (i) by a 
bound-neutron wavefunction, 

-|«»^" 

and (ii) by a quasibound antineutron reduced wavefunction v U £j(r) — w v ij(r) / 5m, 



-u„i 3 (r) + — — u vlj (r) - U n (r)u utj (r) + B^ u vlj (r) = , (7) 



- (r) + 2 l 2 } vutj (r) - [U n (r) + i W* (r)]v vtj (r) + flg v v{j (r) + u vlj (r) = . (8) 

Operating on Eq. ([8|) similarly to the way in which Eq. (J6j) was derived from Eq. J4}, recalling that B ™^ and u u gj(r) 
are real, and multiplying the result by (5m) 2 in order to make connection with Eq. ([6]), we obtain 



- 2 (5m) 2 J Im (v uij (r) u* ut] (r)) dr = 2 (5m) 2 J W n (r) | v uij (r) f dr , (9) 

so that the disappearance width from the vtj sp state, Eq. |6]), is given by 

_ 2 (,5m) 2 J Wn(r) | ^(r) | 2 rfr _ 2 (<5m) 2 J M^ (r) Im ^(r) rfr 

in terms of the solutions u v tj and of Eqs. ([7]) and ((Sj> , respectively. The averaged disappearance width per neutron 
is then given by 



N 

where n„gj is the appropriate number of neutrons in the sp state v£j, N — ^2 n vij is the number of neutrons in the 
decaying nucleus and summation is over the occupied neutron sp states. Since Yd scales as (5m) 2 , hence inversely 
proportional to T 2 n , it is customary to define a reduced lifetime Tr given by 
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which has dimension of inverse time (s 1 ) . The nuclear disappearance lifetime Td is then given by 

T d =^- d =T R T 2 nfl . (13) 

We solve numerically both Eqs. ([7]) and ([8} for neutron and antineutron sp states, respectively. Eq. ([7|) is identical 
with that used in nuclear bound state problems for occupied sp neutron states. Eq. ([8]) is solved for each one of the 
antineutron sp states with reduced wavefunctions v v ij that are generated by the corresponding occupied sp neutron 
wavefunctions u v ij acting as an inhomogeneous source. The solutions u v ij and v u gj serve as input in the integrals 
Eq. (fl"0|) for the widths T^ij . For completeness we note the precise expression for T v (.j , without neglecting contributions 
of order (5m/ B) 2 : 

Tvi . = 2/W ft (r) \w vtj (r) | 2 dr 
V 3 J (I Uuij{r) | 2 + | w v tj{r) | 2 ) dr ' 
where u v gj and w v tj solve Eqs. |4]) and ([5]). 



B. Nuclear structure issues 

The standard shell-model (SM) description of neutron sp states in nuclei introduces spurious excitations of the center 
of mass degree of freedom. It is important, particularly in as light a nucleus as 16 O, to eliminate this spuriosity and thus 
avoid its unphysical effects on the equations solved for the neutron sp states and on the neutron disappearance widths 
subsequently derived. A general construction of translationally invariant (TI) nuclear wavefunctions and densities, in 
the harmonic-oscillator basis, was given by Navratil l3]. Here we follow the earlier discussion by Millener et al. 
which is specifically suited to the p shell. Solving radial equations in the relative coordinate between a neutron and 
its nuclear core in a p-shell nucleus, the number of neutrons in the s shell and p shell has to be modified from the SM 



values nf™ = 2 and nf™ = (N - 2) to 

^ = 2 -^f> n™ = ^(N-2), (15) 

where for N = 8 appropriate to 16 we have nj} = 1.6 and nfp — 6.4. To reproduce a given value of the mean-square 
(ms) radius of the point-neutron distribution < r 2 >„, the ms radii of the neutron-core Is and lp wavefunctions have 
to satisfy 

1 A-\ 
4 

where it was assumed that the s-hole and p-hole strengths are not fragmented. In practice we used a spin-orbit 
potential to split the p-hole strength according to the observed px/2 ~ P3/2 energy difference. Equation p^|) is to be 
contrasted with the SM version in which center of mass spuriosities are disregarded: 

1 



< r * >n= -(— H 2 ( n S < r * >i« + n Tl < r2 >1p) - (16) 



< r 2 > n = ^ (n™ < r 2 > ls + < r 2 > lp ) . (17) 



Numerical solution 



The real wavefunction u„tj for the bound neutron is obtained by solving numerically Eq. ([7]) using a standard 
method. Here we describe briefly the method used to solve the inhomogenous equation ([5]) for the quasibound 
antineutron reduced wavefunction in the potential taken from fits to antiprotonic atom data [10(. 
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The inhomogeneous radial equation (jSJ) is integrated numerically from r — outward using the Numerov method, 
requiring the usual regular boundary conditon v v ij ~ r f+1 at the origin. In parallel, the corresponding homogeneous 
radial equation 

" + ^-f^ffiW - + i Wn(r)]«g(r) + flfg «<»(»■) = • < 18 ) 

obtained from Eq. ((SJ) by omitting the last term, is also integrated using the regular boundary condition at r = 0. 
Integration is carried out to a matching radius R where the nuclear n potentials may safely be neglected. Both 
integrations lead to exponentially increasing functions toward R, as expected. The most general, regular at r = 
solution of Eq. ((SJ) is given by the linear combination 

v uij = <Wj v vlj + v »tj > ( 19 ) 

where the (complex) constant a u tj is chosen such that u^v is regular also at infinity. 

We note that outside the matching radius R the homogeneous equation (fT5|) is satisfied by the neutron bound-state 
wavefunction u v tj . The most general, regular at r — > oo solution of Eq. © is then given by the linear combination 



v 



> 



Kij u vij + v vi j , (20) 



where v u ej is a special, regular at r — > oo solution of the inhomogeneous equation (jHJ), and b u ij is an arbitrary 
(complex) constant. 

The constants a v ij and b„ij above are determined by requiring that the 'inside' and 'outside' solutions, and 
v ^ej respectively, as well as their first derivatives, agree with each other at the matching point r = R. We note that 
for the purpose of evaluating the disappearance widths T„ij, only the knowledge of the constants a U £j is required. A 
straightforward algebra gives 

KiiWutjiR) - u vtj {R)v' uij {R)] - [u' uij {R)v vlj {R) - u vlj {R)v' vlj {R)] , 01 . 

a "t3 = 7Z. • 

< ej (R)v { °i(R) - u v , 3 {R)v^ (R) 

Since both u v ij and v v ij fall off exponentially, whereas v v ij increases exponentially at R, the terms in the first square 
bracket in the numerator of Eq. (|2ip are of order exp(— 2K v ijR) with respect to the terms in the second square bracket, 
and may be safely neglected. Here K u (j — (2/x B^.) 1 / 2 /%, and exp(— 1n v njR) < for a wide range of realistic 

nuclear applications. Hence, to such high accuracy 

<ij (R)vut, (R) - u„ ej (R)v' vej (R) 

<^R)v%(R)~u„ l3 (R)v%(R) 

and no specific knowledge of the special solutions v u ij is required. [Special solutions v u ?j in terms of neutron wave- 
functions u v gj are given in the Appendix.] A further simplification of Eq. (|22[) occurs by noting that v v ij in the 
numerator and in the denominator, for any given t value, share the same asymptotic behavior at R, which leads 
to the final result 

v vij\R) 

These expressions are useful only when their dependence on the matching radius R is negligible. In practice we used 
Eq. J2"2"j) with radial integration steps of 0.04 fm and R between 10 and 13 fm. The coefficients a v ij which determine 
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FIG. 1: Antineutron lpi/2 reduced radial wavefunction v in O. 
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FIG. 2: Integrands in the numerators for the two versions of the neutron disappearance width Eq. (|10|l . 

the required decaying h radial wavefunctions in the region r < R were found to vary by less than 10 -5 of their values, 
with the resulting widths stable to better than 10~ 4 . Both forms of Eq. (fTU|) gave identical results to this order. The 
same accuracy is also obtained using Eq. (|23[) . For comparison with previous calculations, we calculated the widths 
for 16 O, using the input paprameters given by Dover et al. [8|. The results agreed with those listed in Table I of that 
reference to within 1%. 

Figure [T] shows, as an example, the reduced antineutron lpx/2 wavefunction v in 16 0. Note that the longer tail of the 
real part is of no consequence for the disappearance width because in the | v v ij{r) | 2 version of the integral, Eq. (|10p . 
the imaginary potential W is of a shorter range. The other version of the width contains only the imaginary part 
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TABLE I: Calculated reduced neutron disappearance widths T^ij /(8m) 2 (in MeV 1 ). 



state vij 




SM 
n vlj 




ls O(SM) 


16 0(TI) 


56 Fe(SM) 


lSl/2 




2 


1.6000 


0.0207 


0.0220 


0.0144 


1P3/2 




4 


4.2667 


0.0296 


0.0296 


0.0177 


lPl/2 




2 


2.1333 


0.0321 


0.0343 


0.0177 


14/2 




6 


— 


— 


— 


0.0217 


2si/ 2 




2 


— 


— 


— 


0.0238 


ldz/2 




4 








0.0219 


1 J7/2 




B 
o 








u.uzuo 


2p 3 /2 




2 








0.0343 


r d /(5m) 2 


-average [Eq. (JTTJ] 






0.0280 


0.0294 


0.0228 


T d /(5m) 2 


-closure [Eq. ([2j)] 






0.0271 


0.0265 


0.0220 


T R [Eq. { 


H (s- 1 ) 






0.543 xlO 23 


0.517 xlO 23 


0.666 xlO 23 



of the antineutron wavefunction, which is of a shorter range as is seen in the figure. Figure [2] shows the integrands 
in the numerator for the two versions of the width, for lsi/2 and lpi/2 neutrons in 16 0. It is remarkable that the 
relative differences between the integrals are less than 10~ 5 , although the two integrands are not identical both near 
the maxima and in the tail region. 

III. RESULTS 
A. Neutron disappearance widths 

Neutron disappearance widths for the various subshclls vij were calculated for 16 and 56 Fe with the two forms 
given by Eq. (| lOj) . The bound neutron wavefunctions u u ij were calculated in a Woods-Saxon potential — /[l + 
exp((r — i? 1 / 2 )/a)] whose depth was adjusted in each nucleus to fit the experimental separation energy of the 
least-bound neutron. A spin-orbit term was added to reproduce the observed px/2 — P3/2 splitting. The half-density 
radius parameters R\/2 — Tq(A — l) 1 / 3 were adjusted such that the root-mean-square (rms) radius of the whole 
neutron distributions was 2.57 fm for 16 and 3.71 fm for 56 Fe. For 16 it corresponds to the known rms radius for 
the point-proton distribution and for 56 Fe it is 0.09 fm larger than the known value for the point-proton distribution 
[l9l | . The diffusivity parameter a was fixed at a = 0.60 fm for 16 and a — 0.55 fm for 56 Fe. Using Eq. (| 1 T|) for 16 and 
its straightforward extension for 56 Fe, values of tq = 1.325 fm and ro = 1.304 fm were found for the SM calculations 
in 16 and 56 Fe, respectively, whereas using Eq. (TlT))) a value of ro = 1.442 fm was found for the TI calculation in 16 0. 
The depths of the SM potential were 53.8 MeV for 16 and 51.1 MeV for 56 Fe, and the depth of the TI potential in 
16 was 48.8 MeV. For the antineutrons we used the most recent phenomenological isoscalar potential obtained from 
large-scale ('global') fits to 90 data points of strong-interaction shifts and widths in antiprotonic atoms across the 
periodic table [lfj. The effective amplitude for that potential is bo = 1.3 + i 1.9 fm, used with a finite-range Gaussian 
folded with a range parameter 0.9 fm into the nuclear matter density; see Ref. [3] for more details. These potentials 
are shown in Fig. [3] and Fig. [4] for n— 15 and n— 55 Fe, respectively. 
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FIG. 3: The antineutron optical potential in le O. 




FIG. 4: The antineutron optical potential in Fe. 



The neutron disappearance widths as given by Eq. (flO|) scale with (8m) 2 . Therefore in TableUwe list the calculated 
reduced neutron disappearance widths T v £j/(Sm) 2 for the various sub-shells in 16 and in 56 Fe. These calculated 
widths increase by as much as 50% in 16 and by over 100% in 56 Fe going outward from the inner lsx/a sp neutron 
state to the least-bound sp neutron state, in agreement with the trend found in earlier calculations |8|. The values of 
Td/ (6m) 2 calculated using the closure expression Eq. ([2]) give excellent approximation to the exact weighted averages, 
provided values of the nuclear potentials U n , U n , W% at the nuclear surface, here approximated by half of the 
corresponding values at the center of the nucleus, are adopted. We comment that for the strongly absorptive n 
potential used by us, the effect of the real potentials U n and U n is secondary to that of the imaginary potential W n . 
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TABLE II: n scattering lengths a (in fm) from fits to p atomic data. 



Nucleus 


FR \W} 


ZR [10] 


Eq. iPJl 


S [22] 


D [22] 


12 C 


3.26 - i 0.97 


3.21 -i 0.84 


(3.34 ±0.07) -i (1.00 ±0.04) 


3.16 -i 0.87 


3.11 - i 1.10 


58 Ni 


5.41 - i 1.12 


5.43 - i 1.14 


(5.44 ±0.11) -i (1.00 ±0.04) 







A more elaborate averaging with the nuclear density should give similar results for strongly absorptive n potentials, 
as found by Dover et al. [8J]. As for using the TI or the SM schemes for 16 O, the difference between the averaged 
widths, or between the reduced lifetimes, amounts to merely 5% which affects the limit placed on T n „ to only 2.5%. 
Our calculated values of the reduced lifetimes Tr are smaller than those calculated in Ref. [8] . This difference reflects 
partly the difference in the n potentials used and partly the more precise treatment of the nuclear geometry in our 
calculations. The half-density radius parameter of the WS potential used by Dover et al. in the SM calculation for 
16 was taken as R\/2 = 2.545 fm, whereas it assumes the value R\/2 = 3.268 fm in the present SM calculation for 
16 0. The values of R\/2 used for 56 Fe are very close to each other in these works. Finally, for the recently reported 
SK result, 7^(0) > 1.77 x 10 32 yr [jjj in oxygen, our calculated value for the reduced neutron disappearance width, 
using the TI version, implies a lower limit of r nfi > 3.3 x 10 s s. 



B. Antineutron scattering lengths 



For ultracold neutron experiments searching for nn oscillations, the knowledge of their scattering lengths in a given 
material is essential for constructing the relevant Fermi pseudopotentials. The neutron scattering lengths are known 
from thermal neutron reactions. Here we discuss the extraction of antineutron scattering lengths from the p optical 
potentials determined from a comprehensive analysis of p atomic data. We note that for the best-fit p potentials, the 
isovector component came out negligible 



were discussed extensively in Ref. 
form 



lfj, so the n and p potentials are identical. Antinucleon scattering lengths 



2fJ . It was found there that a simple global parameterization was possible, in the 



Re a = (1.54 ± 0.03)^' 



0.311±0.005 



fm , Im a = -1.00 ± 0.04 fm 



(24) 



for A > 10. The approximate A 1 / 3 dependence of the real part, and the constancy of imaginary part, are to be 
expected on the basis of a simple model based on a strongly absorptive square well potential, although the actual 
magnitude of Im a is considerably larger than expected for a sharp-edge potential, resulting mainly from the diffuscness 
of the potential 



21|. 



In our latest work [lOj the data base has been extended to include the numerous CERN PS209 Collaboration p 
atomic data. The values of a due to the isoscalar N potentials fitted to this extended set of data are listed for 12 C and 



22]. Here, 



58 Ni in Table ITT1 together with values obtained using Eq. (|24|) . and also from the earlier work by Wong et al 
FR and ZR stand for values of a derived from best global-fit finite-range and zero-range potentials, respectively, with 
FR giving the lowest x 2 [101 ■ The notations S (shallow) and D (deep) stand for values of a calculated in Ref. [15] from 



potentials with a relatively weak absorptivity W or a strong one, respectively, derived from limited fits in Ref. 22 1. 
The resulting scattering lengths a exhibit a remarkable independence of the model used, provided it fits the p atomic 
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data. This stability follows from the strong absorptivity of the p potential which suppjresses the associated Is atomic 
radial wavefunction in the nuclear interior where the main model dependence arises 



3- 



IV. SUMMARY 



We have reported results of precise calculations of the reduced lifetimes of representative nuclei, 16 O and 56 Fe, 
against neutron-antineutron oscillations, thus providing revised and updated lower limits on the free-space nn, oscil- 
lation time T nfl . The best lower limit is now provided by the very recent SK measurement [l^ | in 16 O which yields 
according to our calculation a limit of T nfl > 3.3 x 10 8 s. We have used the latest (isoscalar) antinucleon potentials 
derived from the analysis of a large-scale set of p atomic data near threshold Having solved accurately the sp 
equations for neutrons and (coupled) antineutrons, it became possible to test the usefulness of rough approximations 
such as Eq. ([2]) in terms of mean nuclear potentials for neutrons and antineutrons. We found that using surface values 
for these mean potentials, taken as half the corresponding values in the center of the nucleus, provided an excellent 
approximation to the exact calculation. This points out to a considerably mild model dependence of the calculated 
reduced widths that are sensitive foremost to n, potentials at the nuclear surface where their determination from p 
atomic data involves only little extrapolation. 

An educated estimate of the theoretical uncertainty involved in the derivation of the lower bound deduced in the 
present work on the nn oscillation time r nfi can be made as follows. 

• For a given nucleus like 16 O and a given n-nuclear potential, but considering alternative ways of treating the 
nuclear size, the calculated nuclear widths vary by 5% (columns 4 and 5 of Table [T| , so the uncertainty in the 
derived r nfl is about 2.5%. 

• The uncertainty arising from using different nuclei (columns 4 and 6 of Table [TJ) comes mostly from the uncer- 
tainty in the strength of the absorptive (imaginary) n-nuclear potential. That shows about 20% uncertainty for 
the averaged disappearance width T ( i, and hence 10% for T n n. 

Thus, the overall theoretical uncertainty involved in the present one-nucleon nn oscillation calculations is about 
10% — 15%. It should be viewed as a model-dependence uncertainty that is considerably lower than the 50% — 100% 
uncertainty range evident in many of the calculations from the 1980s and 1990s, e.g. Ref. [8j, before the information 



lCj. Other past calculations 



flfl 



which 



from p atoms became as abundant and precise as it is available to date 
avoided using h phenomcnological optical potentials faced a tougher task of having to renormalize the nN strong 
interaction within the nucleus, a formidable job that was bypassed by Dover et al. [8[ and in the present work using 
a well-constrained phenomenological ri-nuclear potential. 

Another source of uncertainty involves two-nucleon processes which inside the nucleus might compete with the 
leading one-nucleon process considered here. In their 1985 paper, making contact with beta-decay and EMC calcula- 
tions, Dover et al. [8] estimated these additional modes of neutron disappearance to be about 15% — 30%, and largely 
incoherent with the one-nucleon mode. This provides a systematical uncertainty which may be used to increase the 
stated lower bound on r™ . 
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Appendix: asymptotic n radial wavefunctions 



Here we record special, regular at r — * oo solutions of the n inhomogeneous radial Eq. (jHJ) in terms of similar 
solutions of the n homogeneous radial Eq. ([7]). At r > R, where the n and n nuclear potentials are negligible, these 
equations are written in standard form as 

£fg _|_ 1) 1 

- v utj {p) + j — V ^j(p) + v v tj{p) + —nxu„ ej (p) = , (25) 

9 B ulj 
£(£ _|_ 1) 

- u vVj{p) H ~2 — u *v(p) + u vlj{p) = , (26) 

for n and n, respectively, where p = K u tjT is dimensionless. The n bound-state, regular at p — > oo solutions are given 
by 

u vli {9) = A uej V e (p)eM-p) = A vl] {-1) 1 p e+1 (-^-) e ^tfl ; (27) 

p dp p 

where A u gj are normalization constants ensuring that asymptotically Vi(p — > oo) — > 1. Ve(p) are polynomials in 1/ p, 
related to the outgoing spherical Hankel functions 24|, satisfying the differential equation 



V'l{p) - 2V' e (p) - f ^±llp e (p) = . (28) 
P 

The lowest-order Vgs relevant to the present work are 

V (p) = l, V 1 (p) = (l + -), V 2 { P ) = (l + - + \) , ^3(^)^(1 + ^ + ^ + ^). (29) 

p p p p p p 

A useful recursion relation satisfied by the Ves is 

Ve(p) = (l + -)Ve- 1 (p)-V' i _ 1 (p) (P_i(p) = l), (30) 
P 

easily derived from the explicit form of Vg given in Eq. (|27[) . 

It can be shown that regular at p — > oo solutions of the n inhomogeneous radial Eq. (|25[) are given by 



A p ■ 

Vvtj{p) = t^t P 2Kf)exp(-p) , (31) 

where Qe(p) = Vt—i{p)- For a proof, one forms the inhomogeneous second-order differential equation satisfied by the 
Qis, making use of Eqs. (f2"8")) and ([3TJ]) . The latter equation allows us then to construct the Qgs recursively: 

Qe +1 (p) = (l + - p )Qdp)-Q e (p) (Qo(p) = l). (32) 
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